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Abstract. We investigate the origin of color gradients 
in cluster early-type galaxies to probe whether pure age 
or pure metallicity gradients can explain the observed 
data in local and distant (z w 0.4) samples. We mea- 
sure the surface brightness profiles of the 20 brightest 
early-type galaxies of CL0949-I-44 (hereafter CL0949) at 
redshift z=0.35-0.38 from HST WF2 frames taken in the 
filters F555W, F675W, F814W. We determine the color 
profiles {V — R){r), {V — I){r), and {R — I){r) as a func- 
tion of the radial distance r in arcsec, and fit logarithmic 
gradients in the range —0.2 to 0.1 mag per decade. These 
values are similar to what is found locally for the colors 
{U — B), {U — V), {B — V) which approximately match 
the {V -R),{V -I),{R- I) at redshift w 0.4. We anal- 
yse the results with up to date stellar population mod- 
els. We find that passive evolution of metallicity gradients 
(« 0.2 dex per radial decade) provides a consistent ex- 
planation of the local and distant galaxies' data. Invoking 
pure age gradients (with fixed metallicity) to explain local 
color gradients produces too steep gradients at redshifts 
z w 0.4. Pure age gradients are consistent with the data 
only if large present day ages > 15 Gyr are assumed for 
the galaxy centers. 

Key words: galaxies: clusters: general - galaxies: ellip- 
tical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: fundamental parameters 



1. Introduction 

The efforts to understand the mechanisms of formation 
and evolution of elliptical galaxies have dramatically in- 
creased in the last years. On the theoretical side, the clas- 
sical picture of the monolithic collapse a la Larson ( 1974 ) 
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has been increasingly questioned by scenarios of hierar- 
chical structure formation, where ellipticals are formed 



by mergers of spirals (White & Rees 1978). While in the 



first case the stars of the galaxies are supposed to form 
in an early episode of violent and rapid formation, the 
semi-analytical models of galaxy formation (Kauffmann 



1996 , Baugh et al. 1996 ) predict more extended star for- 
mation histories, with substantial fractions of the stellar 
population formed at relatively small redshifts. From an 
observational point of view, evidence is growing that the 
formation of a fair fraction of the stellar component of 
bright cluster Es must have taken place at high redshift 
(z > 2) and cannot have lasted for more than 1 Gyr. 
This conclusion is supported at small redshifts by the ex- 
istence of a number of correlations between the global pa- 
rameters of elliptical galaxies, all with small scatter: the 
Color-Magnitude (CM, Bo wer et al. |1992| ), the Fundamen- 
tal Plane (J0rgensen et al. 199(:), the Mg-cr (Colless et al. 



1999) relations. These hold to redshifts up to 1 (Stan 
Bender et 



1998 



al. 



1998, Bender, Ziegler 



1998| , Kelson et al. 
Bruzual 



1997 



1997, Ellis et al. 1997), consistent with 



1996, 



ford et al 
van Dokkum et al 
Ziegler & Bender 
passive evolution and high formation- ages. In addition, the 
Mg over Fe overabundance of the stellar populations of lo- 
cal cluster Es (Worthey, Faber & Gonzalez 
et al. 



1992 



Mehlert 



2GGGb), requires short star formation times and/or 
a flat IMF (Thomas, Greggio, Bender |l99|). Field ellipti- 



cals m ight have larger spreads in formation ages (Gonzales 
1993 ), but comparisons with cluster Es (Bernardi et al. 
1998) and the first trends observed for the evolution of 



the FP with redshift (Treu et al. |1999| ) do not allow much 
freedom. 

The study of the variations of the stellar populations 
inside galaxies, the so-called gradients, offers an addi- 
tional tool to constrain the models of galaxy formation. 
Classical monolithic collapse models with Salpeter IMF 



(Carlberg 1984) predict strong metallicity gradients. The 



merger trees generated by hierarchical models of galaxy 



formation (Lacey & Cole 1993) are expected to dilute 



gradients originally present in the merging galaxies and 
therefore give end-products with milder gradients (White 
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198C), although detailed quantitative predictions are still 
lacking. Color gradients are known to exist in ellipticals 
since the pioneering work of de Vaucouleurs ( |1961| ) . The 
largest modern, CCD based, dataset has been collected 
by Fra nx et al. ( |1989| ), Peletier et al. ( |1990| ), Goudfrooij 



et al. (1994), who measured surface brightness profiles in 
the optical U, B, V, R, I bands for a large set of lo- 
cal ellipticals. These radial changes in colors have been 
tradi tionally interpreted as changes in metallicity (Faber 



1977 ), arguing that line i ndex gradients measu red in local 



field (Davieset al. 1993 , CaroUo et al. |l993, Kobayash i 
& Arimoto |l999[ ) and cluster (Mehlert et al. |l998| , ^000a|) 
galaxies support this view. A w 0.2 dex change in metal- 
licity per radial decade suffices to explain the photomet- 
ric and spectroscopic data (see, for example, Peletier et 
al. 1990). However, due to the age- metallicity degener- 



acy o f the spectra of simple stellar populations (Worthey 
1994) this conclusion can be questioned. De Jong ( 1996|) 



favors radial variations in age to explain the color gra- 
dients observed in his sample of local spiral galaxies. If 
ellipticals form from merging of spirals, age might be the 
driver of color gradients in these objects too. The com- 
bined analysis of more metal- sensitive (i.e., Mg2) and more 
age-sensitive (i.e., H(3) line index profiles points to a com- 
bination of both age and metallicity variations (Gonzales 
1994 Mehlert et al. |2000bD Fina lly, Goudfrooij & de Jong 
(1995D and Wise & Silva (1996) stress the importance of 



dust on the interpretation of color gradients and conclude 
that it might add an additional source of degeneracy to 
the problem. 

With the resolution power of HST it is now possible to 
measure the photometric properties of distant galaxies and 
therefore implement an additional constraint on the stel- 
lar population models, i.e. the time evolution, that helps 
break ing the degeneracies discussed above. Abraham et al. 
(1999) study color distributions of intermediate redshift 
galaxies of the Hubble Deep Field to constrain the star 
formation history of the Hubble Sequence of field galax- 
ies. Here we determine the color gradients in the V, R, I 
Johnson-Cousins filters of a sample of cluster early-type 
galaxies at redshift z w 0.4 observed with HST, and com- 
pare them to the U, B and V gradients observed in local 
Es. By means of new, up to date models of Simple Stel- 
lar Populations (SSP), we investigate whether pure age or 
pure metallicity gradients can explain the observed colors 
and gradients and their evolution with redshift. The struc- 
ture of the paper is the following. In Sect. we present the 
new HST observations of CL0949, the data reduction and 
the derived photometric profiles. In Sect. 0we describe the 
reference sample of local cluster galaxies. In Sect. H we de- 
scribe the stellar population models and their use to study 
the evolution of color gradients. Ages and metallicities de- 
rived according to the different modeling assumptions are 
presented in Sect. H Results are discussed in Sect. O and 
conclusions are drawn in Sect. R. 



2. HST Observations of CL0949 

The central core of CL0949 (a superposition of two clus- 
ters at redshifts 0.35 and 0.38, Dressier & Gunn 1992 ) 
was observed with HST and WFPC2 in December 1998 
(PID 6478). Nine images with each of the filters F555W, 
F675W, and F814W were obtained totaling exposure 
times of 7200 sec, 7200 sec, and 7000 sec respectively. 
The images were taken with sub-pixel shifts to improve 
the resolution and sampling. The pipeline reduced frames 
were further processed under IRAF to eliminate cosmic 
rays, and combined after resampling to half the origi- 
nal scale, following the algorithm described in Seitz et 
al. (1998). The subsequent analysis was performed under 
MIDAS. Particular care was needed to treat the images 
taken with the filter F555W, which suffered from an ele- 
vated background level throughout the field of view, due 
to Earth reflection by the optical telescope assembly and 
spiders, and the related presence of a dark diagonal X 
structure. In order to determine the background accu- 
rately, the Sextractor program (Berlin & Arnoux 1996 ) 
was used. We detected and subtracted sources iteratively 
and constructed an average source-free sky frame. The 
resulting source-free, sky-subtracted images taken with 
the F675W and F814W filters have rms residuals of 0.3- 
0.5%, the ones with the F555W filter up to 1%, which 
set the final precision of our sky subtraction. In addi- 
tion, we constructed background-subtracted images using 
the first-iteration background determined by Sextractor, 
which takes into account the diffuse light contribution of 
extended objects. These frames were then used to perform 
the surface photometry analysis for all galaxies except for 
the cD on the WF2 frame (#201 of Fig. |l]). In this case we 
preferred the sky-subtracted images to be able to measure 
the extended halo of the object. 

The isophote shape analysis was performed follow- 
ing Bender and MoUenhoff (1987). The procedure was 
improved in two aspects. 1) Mask files were produced 
automatically from the "segmentation frames" produced 
by Sextractor. 2) The surface brightness and integrated 
magnitude profiles following elliptical isophotes were com- 
puted from the averages of the not masked pixels between 
the isophotes. Statistical errors were also derived from 
the measured rms. Profiles of surface brightness, elliptic- 
ity, position angle and higher order terms of the Fourier 
decomposition were derived for the 70 brightest objects 
detected by Sextractor on the F814W frames. The infor- 
mation was used in combination with the visual inspec- 
tion of the images to classify the galaxies morphologi- 
cally. Twenty early- type galaxies were found. Their iden- 
tification numbers can be read in Fig. H. In the follow- 
ing we shall consider only these objects. The data for the 
late-type galaxies will be presented elsewhere. The cross- 



identification with Dressier & Gunn (1992) is given in 
Table ^ where the available spectroscopic information is 
also summarized. Inspection of the proflles of the isophote 
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Fig. 1. The F814W image of the cluster CL0949. The la- 
bels indicate the 20 brightest early-type galaxies detected. 

shape parameters shows smooth variations as a function 
of the radial distance, and absence of strong odd-term co- 
efficients, which excludes the presence of large amounts of 
patchy dust. This is confirmed by the examination of the 
two-dimensional (V-R), (V-I), (R-I) color maps. 

Circularly averaged surface brightness profiles were 



also measured following Saglia et al. (1997a). The pho- 
tometric calibration on the V, R, I Johnson-Cousins 
bands was p erform ed using the "synthetic syste m" (H olz- 
mann et al. |l995|) as in Fig. 1 of Ziegler et al. (|l999|) . A 
first estimate of the zero-points was obtained using the 
{V — I) and {V — R) colors typica l of z — 0.4 early- type 
galaxies, as in Ziegler et al. ( 1999| ). The final zero-points 
were derived recomputing the calibration with the so- 
measured (V — I) and {V — R) colors. Half-luminosity radii 
iie(555), i?e(675), i?e(814) and total magnitudes tot(555), 
tot(675), 771^(814) were derived fitting the circularly av- 
eraged profiles using the algorithm described and exten- 
sively tested in Sagha et al. ( 1997b ). The algorithm fits 
a PSF (computed using the Tinytim program) broadened 
r^''^ and an exponential component simultaneously and 
separately to the circularly averaged surface brightness 
profiles. As a result, disk-to-bulge ratios and the scale- 
lengths of the bulge and the disk components are also de- 
termined. The quality of the fits were explored by Monte 



Carlo simulations in Saglia et al. (1997b), taking into ac- 
count the possible influence of sky subtraction errors, the 
signal-to-noise ratio, the radial extent of the profiles and 
the x^ quality of the fit. Accordingly, the half- luminosity 
radii are expected to be accurate to 25%, the total magni- 
tudes to 0.15 mag. Two fits are more uncertain than these 



figures. The values of i?e(814) and i?e(675) of galaxy 201 
have been derived restricting the fit to the radial range of 
the F555W profile. Fitting the whole F814W profile, more 
extended than the F555W and F675W ones and shallower 
than the extrapolated profiles expected from the fits in 
the V and R bands, one would get a value of i?e(814) 
twice as large. Fitting the whole F675W profile produces 
a value of i?e(675) 30% smaller. Finally, the F555W profile 
of the galaxy 408 is not extended enough to allow a reli- 
able determination of i?e(555). The best fit the these data 
involves a large (50%) extrapolation and gives a value of 
i?e(555) twice as large as i?e(675) and i?e(814). 

Fig. g shows the resulting fits to the F814W filter pro- 
files and the profiles in the F675W and F555W bands. 

Half-luminosity radii and total magnitudes are given 
in Table |l|. The scales (ReB and h in arcsec) and the mag- 
nitudes {iTiB and itid) of the bulge and disk components 
measured with the F814W profiles are given in Table g 
together with the D/B ratio, for the objects where a two- 
component model was fit. As described above, the fit to 
the galaxy 201 is highly uncertain. If the whole radial ex- 
tent of the F814W profile is fit, an 8 times larger value of 
the bulge scalelength is needed. 



Following Burstein & Heiles (1984), no correction for 
galactic absorp tion was necessary. Note, however, that 
Schlegel et al. ( |1998D give E{B -V) = 0.019. The use of 
this correction gives Ay = 0.063, Ar = 0.051, Ai = 0.037, 
making the derived intrinsic (V — i?), (F — /), (R — I) col- 
ors bluer by 0.01, 0.026, 0.014 mag respectively. These 
variations do not affect what follows. 

Finally, we determined the color profiles of the galax- 
ies from the calibrated elliptical isophote surface bright- 
ness profiles, estimating statistical and systematic errors. 
The latter were computed considering the extreme cases 
of possible sky subtraction errors, ±0.5% for the F675W 
and F814W filters, and ±1% for the F555W (see above). 
For example, we computed the two {V — T) profiles cor- 
responding to -|-1% sky error on V and -0.5% sky error 
on I, and -1% sky error on V and +0.5% sky error on I. 
The error budget on the colors is dominated by statistics 
in the central regions of the galaxies, and by systematics 
in the outer parts. We fitted the logarithmic slope gradi- 
ent and zero point, weighting the data according to the 
statistical errors. In the fit we consider only points within 
0.1 (i?mi„)and 3.2 {Rmax) arcsec and with statistical and 
systematic errors less than 0.1 mag. Systematic errors on 
the gradients and zero-points due to sky subtraction were 
estimated determining the two quantities for the two ex- 
treme cases of possible sky errors discussed above. The 
median value of Rmin/ Re is 0.2, of Rmax/ Re 1-6 for the 
{V - R) and iy - I) colors, and 2.8 for {R - /). We ex- 
plored the effect of the radial range by repeating the fits 
on a fractional (0.2 < R/Re < 1.5) or on a fixed radial 
range (0.1 < i? < 1 arcsec). In both cases we obtained 
similar gradients and zeropoints within the statistical er- 
rors, if the systematic errors in this radial range due to 
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Fig. 2. The fits to the circularly averaged surface brightness profiles (taken with the F814W filter and calibrated in 
the Cousins I band) of the 20 brightest early- type galaxies of CL0949. Two plots arc shown for each galaxy. The upper 
panel shows the calibrated surface brightness profile (small dots) together with the fitted bulge and/or disk profiles 
(full lines) as a function of the 1/4 power of the distance r from the center in arcsec. The horizontal dashed line shows 
the position corresponding to 1 per cent of the sky level. We also show the circularly averaged surface brightness 
profiles taken with the F675W and F555W filters, calibrated in the Johnson-Cousins V and R bands. For clarity they 
have been shifted by 0.5 and 1 mag respectively. In each panel the galaxy name and the type of best fit are given (B 
for bulge only, D for disk only, BD for bulge-f disk fit), together with the value of Re in arcsec and the total magnitude 
calibrated in the Cousins I band. The Re location is also marked as a vertical line from the top xaxis in each panel. 
Below each surface brightness profile we show the residuals of the fits, A^{r) — jiobsir) — M/ii(^) (dots, scale -0.2 to 
0.2 mag arcsec"^) and the integrated magnitude difference Ama(7(r) — magobs{< f) — magfit{< r) (full lines, scale 
-0.1 to 0.1 mag). 



Table 1. Results of the photometric fits. Column 1 gives the galaxy name (see Fig. |l|), Col. 2 the cross identification 
with Dressier & Gunn (1992), Col. 3 the redshift from Dressier & Gunn (1992), Col. 4 the morphological type, Cols. 
5-10 give the half-luminosity radii Re in arcsec and total magnitudes m^ in the F555W, F675W, and F814W filters 
respectively. Magnitudes are calibrated to the V, R, I Johnson-Cousins system (see text). Col. 11 gives the central 
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Name 


D&G92 


z 


Type 


i?e(555) mT(555) 


i?e(675) 


mT(675) 


i?e(814) 


mT(814) 


Comment 


201 


118 


.3815 


cD 


1.46 19.44 


1.90" 


18.06" 


2.05" 


17.23" 


cr = 325 ± 50 


202 






E 


0.64 20.54 


0.59 


19.53 


0.55 


18.78 




203 


102 


.3766 


SBO 


0.35 20.52 


0.37 


19.58 


0.36 


18.87 


E+A 


204 


130 


.3402 


E 


0.63 20.77 


0.58 


19.85 


0.55 


19.09 




205 






E/SO 


0.52 20.91 


0.49 


19.86 


0.48 


19.08 




206 






E 


0.72 20.84 


0.58 


19.91 


0.54 


19.17 




208 


123 


.3779 


SO 


0.45 21.51 


0.47 


20.4 


0.50 


19.58 




209 


110 


.3744 


S/SO 


0.67 21.3 


0.66 


20.34 


0.62 


19.74 




210 






E/SO 


0.44 21.45 


0.34 


20.51 


0.41 


19.69 




211 


108 


.3799 


E/SO 


0.20 21.87 


0.18 


20.79 


0.20 


20.03 




302 


163 


.3778 


E/SO 


0.65 20.42 


0.58 


19.47 


0.57 


18.70 




303 


119 


.3484 


E/SO 


0.43 21 


0.49 


19.88 


0.44 


19.17 




304 






SO 


0.56 21.1 


0.47 


20.11 


0.48 


19.31 




312 






E/SO 


0.15 22.5 


0.16 


21.39 


0.15 


20.67 




401 


221 


.3772 


E 


1.07 20.05 


1.05 


18.87 


1.07 


18.06 


cr = 230 ± 25 


402 


193 


.3767 


Sa/SO 


1.12 20 


1.10 


18.93 


1.12 


18.14 


cr = 230 ± 25 


403 


217 


.3784 


E 


0.69 20.42 


0.71 


19.32 


0.69 


18.57 


a = 125 ±15 


404 


187 


.3778 


E 


0.81 20.63 


0.77 


19.61 


0.76 


18.87 




405 






SBO 


0.37 21.79 


0.29 


20.83 


0.29 


20.05 




408 






SO 


_2 _2 


1.03 


20.7 


0.85 


19.78 


Background? 



1: uncertain 

2: profile too short to allow a reliable fit 



Table 2. The parameters of the bulge-fdisk fits as mea- 
sured with the F814W profiles. Col. 1 gives the galaxy 
name (see Fig. |l|), Col. 2 the half-light radius of the bulge 
component R^b in arcsec, Col. 3 the magnitude of the 
bulge component tub , Col. 4 the scalelength h of the disk 
component in arcsec, Col. 5 the magnitude of the disk 
component m^i. Col. 6 the disk-to-bulge ratio D/B. Mag- 
nitudes are calibrated to the I Cousins filter (see text). 



Name i^^s (814) 


ms (814) 


ft(814) 


mi3(814) 


D/B 


201" 


0.80 


17.81 


3.08 


18.18 


0.71 


203 


0.44 


19.32 


0.18 


20.05 


0.51 


208 


0.16 


20.5 


0.48 


20.18 


1.34 


209 


0.26 


21.42 


0.42 


20 


3.69 


211 


0.18 


20.4 


0.14 


21.36 


0.42 


302 


0.49 


18.89 


0.54 


20.68 


0.19 


303 


0.60 


19.41 


0.13 


20.92 


0.25 


304 


0.38 


19.81 


0.36 


20.41 


0.58 


312 


0.12 


21.13 


0.11 


21.83 


0.53 


402 


1.02 


18.4 


0.80 


19.82 


0.27 


403 


0.67 


18.67 


0.45 


21.21 


0.097 


405 


0.25 


20.48 


0.19 


21.26 


0.49 


408 


0.82 


20.09 


0.54 


21.3 


0.33 


1: uncertain 


fit 











sky subtraction are not too large (« 0.1 mag). We verified 
that the differences between the HST PSFs in the different 
filters do not affect the derived gradients within the errors. 
Tests on i?^/"* profiles of R^ — 0.5 arcsec convolved with 
the PSF for the fihers F555W and F675W show that the 
systematic errors produced on the logarithmic gradients 
computed following our procedure are smaller than the 
systematics due to the sky subtraction errors. For galaxy 
312, the smallest of the sample {Re ~ 0.15"), we derived 
color gradients by first convolving the pairs of images (for 
example, F675W and F555W) with the corresponding in- 
terchanged PSFs (F555W and F675W). The differences in 
the derived logarithmic slopes are smaller than twice the 
statistical error; thus we expect that this effect is negli- 
gible for the other, more extended objects. Figs, p^ to ^ 
show the color profiles with the logarithmic gradient fits. 
Tables 0-0 summarize the results of these fits. The colors 
given there are computed at half Re in the F675W band. 
The galaxy 408 has colors 0.2-0.5 mag redder than the rest 
of the sample. We assume that it is a background object 
and do not consider it in the following analysis. Inspection 
of Figs, y to ^ shows that on the whole the color profiles 
are well described by the logarithmic fits, with mean rms 
residuals in the range 0.02-0.03 mag. However, the fits are 
not optimal in a statistical sense, because their reduced 
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y^ are rather large, in the range 2-4. This stems from the 
small radial scale color variations present in many profiles. 
Similar effects are observed in local galaxies (i.g., Peletier 
et al. [1990D. 



Table 4. The results of the fits to the (V-I) color profiles. 
Columns as in Table 0. 



Table 3. The results of the fits to the (V-R) color profiles. 
Col. 1 gives the galaxy name (see Fig. |^), Col. 2 the fitted 
logarithmic gradient in mag per radial decade, Col. 3 the 
statistical error, Col. 4 the systematic error, Col. 5 the 
fitted color at 0.5 -Re (675), Col. 6 the statistical error. 
Col. 7 the systematic error. 



Name a(V-R) da^t 



das 



Col dColst dCol, 



201 
202 
203 
204 
205 
206 
208 
209 
210 
211 
302 
303 
304 
312 
401 
402 
403 
404 
405 
408 



-0.101 
-0.179 
-0.007 
-0.083 
-0.120 
-0.198 

0.024 
-0.245 
-0.072 
-0.083 

0.026 
-0.067 
-0.155 
-0.058 

0.013 
-0.048 
-0.154 
-0.123 
-0.011 
-0.359 



0.023 
0.019 
0.028 
0.018 
0.035 
0.023 
0.052 
0.044 
0.059 
0.088 
0.033 
0.028 
0.039 
0.066 
0.015 
0.017 
0.011 
0.027 
0.053 
0.063 



0.025 
0.039 
0.035 
0.039 
0.036 
0.043 
0.049 
0.048 
0.050 
0.057 
0.038 
0.043 
0.046 
0.082 
0.034 
0.033 
0.033 
0.040 
0.054 
0.046 



1.178 
1.082 
0.947 
1.019 
1.141 
1.100 
1.139 
0.981 
1.134 
1.155 
1.099 
1.082 
1.144 
1.125 
1.192 
1.115 
1.084 
1.037 
1.131 
1.268 



0.015 
0.016 
0.027 
0.015 
0.031 
0.019 
0.047 
0.033 
0.060 
0.106 
0.026 
0.024 
0.033 
0.083 
0.009 
0.011 
0.009 
0.019 
0.055 
0.046 



0.019 
0.011 
0.002 
0.012 
0.009 
0.014 
0.012 
0.018 
0.005 
0.012 
0.009 
0.010 
0.009 
0.020 
0.017 
0.019 
0.012 
0.016 
0.002 
0.040 



Fig. @ summarizes the properties of the gradient sam- 
ple. No correlation is observed between galaxy size and/or 
galaxy magnitude and gradients. All objects have detected 
negative gradients in at least one color. For at least half of 
the sample the gradients are detected with at least Scr sig- 
nificance. None of the positive gradients is statistically sig- 
nificant. The systematic errors are on average larger than 
the statistical ones. Due to the rather large errors, the 
correlation between the gradients in the different bands is 
not very strong. 

3. The Local Sample 

Since the F555W, F675W, F814W bands at z « 0.4 map 
approximately the [/, i3, V Johnson filters at z = 0, the 
direct comparison of our measured color gradients with 
those derived for local early-type galaxies in the U, B, V 
bands can give a first idea of their evolution. Therefore, 
we consider the s urface brightnes s profi les measured by 
Franx et al. ( |l989D , Peletier et al. ( |l990|) , and Goudfrooij 



Name a(V-I) da^t 



das 



Col dColst dCols 



201 
202 
203 
204 
205 
206 
208 
209 
210 
211 
302 
303 
304 
312 
401 
402 
403 
404 
405 
408 



-0.114 
-0.225 
-0.077 
-0.165 
-0.090 
-0.217 

0.026 
-0.325 
-0.015 
-0.141 
-0.040 
-0.087 
-0.085 

0.032 
-0.008 
-0.077 
-0.207 
-0.182 
-0.037 
-0.616 



0.011 
0.016 
0.041 
0.015 
0.024 
0.031 
0.034 
0.043 
0.039 
0.042 
0.022 
0.018 
0.041 
0.044 
0.010 
0.015 
0.012 
0.031 
0.031 
0.071 



0.028 
0.034 
0.032 
0.037 
0.034 
0.041 
0.046 
0.045 
0.046 
0.052 
0.038 
0.042 
0.043 
0.077 
0.032 
0.032 
0.037 
0.038 
0.051 
0.054 



1.963 
1.864 
1.663 
1.777 
1.896 
1.858 
1.926 
1.688 
1.879 
1.904 
1.862 
1.830 
1.927 
1.797 
1.991 
1.881 
1.835 
1.783 
1.917 
2.241 



0.007 
0.013 
0.039 
0.012 
0.022 
0.026 
0.031 
0.033 
0.041 
0.051 
0.016 
0.016 
0.035 
0.055 
0.007 
0.009 
0.009 
0.022 
0.032 
0.050 



0.021 
0.010 
0.003 
0.012 
0.009 
0.013 
0.011 
0.017 
0.005 
0.010 
0.008 
0.010 
0.008 
0.018 
0.016 
0.018 
0.013 
0.015 
0.002 
0.043 



Table 5. The results of the fits to the (R-I) color profiles. 
Columns as in Table |[ 



Name a(R-I) das 



das 



Col dCols 



dCols 



201 


-0.001 


202 


-0.055 


203 


-0.069 


204 


-0.049 


205 


-0.050 


206 


-0.012 


208 


0.069 


209 


-0.054 


210 


-0.019 


211 


0.149 


302 


-0.024 


303 


-0.059 


304 


-0.002 


312 


0.113 


401 


-0.005 


402 


-0.021 


403 


-0.072 


404 


-0.031 


405 


-0.112 


408 


-0.208 



0.009 
0.016 
0.016 
0.013 
0.028 
0.012 
0.024 
0.015 
0.026 
0.059 
0.024 
0.020 
0.023 
0.058 
0.010 
0.008 
0.010 
0.012 
0.034 
0.022 



0.014 
0.029 
0.027 
0.032 
0.033 
0.029 
0.035 
0.029 
0.034 
0.042 
0.029 
0.034 
0.035 
0.054 
0.023 
0.021 
0.023 
0.031 
0.038 
0.039 



0.788 
0.782 
0.718 
0.757 
0.751 
0.758 
0.790 
0.716 
0.753 
0.673 
0.766 
0.755 
0.779 
0.664 
0.804 
0.769 
0.747 
0.751 
0.808 
0.999 



0.005 
0.012 
0.015 
0.010 
0.023 
0.010 
0.020 
0.011 
0.025 
0.070 
0.018 
0.016 
0.018 
0.072 
0.006 
0.005 
0.007 
0.008 
0.034 
0.013 



0.009 
0.005 
0.000 
0.005 
0.004 
0.006 
0.004 
0.009 
0.000 
0.013 
0.003 
0.003 
0.002 
0.018 
0.009 
0.010 
0.006 
0.009 
0.004 
0.020 



et al. ( 1994 ) . Focusing on cluster galaxies our sample con- 



sists of 21 objects, for which at least one color gradient is 
available. Table pi lists the selected cluster galaxies, their 
membership, the color gradient, the colors at i?e/2 and the 
source of the photometry. For those objects with source of 
photometry "I", the V profile was not available. We de- 
rived it by interpolation of the B and R profiles using the 



I 




-1 -0.5 0.5 



log r(arcsec} 



Fig. 3. The F555W-F675W color profiles (with statistical errors) and the logarithmic gradient fits of the 20 brightest 
early-type galaxies of CL0949. Notice that since the colors of N408 are particularly red, we plot the observed colors 
shifted by -0.2 mag. The shift of 0.2 mag indicated in the figure has to be added to get the real colors of the object. 
The dotted curved lines show the systematic errors on the colors due to sky subtraction. The straight full lines show 
the fitted gradients and their statistical errors. The straight dashed lines show their systematic errors. Relevant radial 
distances are marked on the top x axis of each panel as vertical lines: the location of i?e(675) (solid); the minimum 
and maximum radius used in the fit (dotted); the 0.2_Re — 2i?e range (arrows) used in the modeling. Finally, each panel 
is labeled with the galaxy name, and with the gradient slope (upper row) and zero point (lower row) followed by their 
statistical and systematic error. 



relation {B — V) = 0.66 x {B — R)— 0.03, which was deter- older than 4 Gyr, independently of the metallicities. The 
mined from the SSP models of Maraston (2000) for ages errors on the quantities given in Table ^ are systematic 
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-1 -0.5 0.5 



log r(arcBec) 



Fig. 4. The F555W-F814W color profiles and logarithmic gradient fits of the 20 brightest early-type galaxies of CL0949. 
Lines and labels as in Fig. S. The shift to be apphed to the colors of N408 is shown as in Fig. pi 



and due to the uncertainties in the sky subtraction. The 
statistical errors are (usually) much smaller. Inspection of 
Table |6| shows that our local galaxy sample consists of 16 
objects in {U - B) and {U -V), and 21 in {B -V). In 
the following we indicate with DS the distant sample of 
E galaxies of CL0949 at 2; « 0.4 and with LS the local 
sample of cluster E galaxies from the literature. Figure 
1^ shows the comparison of the cumulative distributions 
of the absolute total B band magnitudes of the two sam- 



ples, having applied 0.4 mag dimming (cf. Ziegler et al. 
|1999| ) to the DS galaxies. The distributions look similar 
(the Kolmogorov-Smirnov probability is 0.28), although 
the LS suffers from incompleteness at magnitudes fainter 
than K, -21. On the whole, the comparison between the 
two samples is fair. 

Fig. compares the cumulative distributions of the 
measured gradients for the LS and the DS samples. The 
distributions appear rather similar, suggesting little evo- 



I 




log r(arcsec) 



Fig. 5. The F675W-F814W color profiles and logarithmic gradient fits of the 20 brightest early-type galaxies of CL0949. 
Lines and labels as in Fig. S. The shift to be apphed to the colors of N408 is shown as in Fig. pi 



lution. In particular, the median values of the {U — B), 
[U — V) and {B — V) gradients of the LS galaxies are 
-0.12, -0.17 and -0.06 mag/dex per decade respectively. 
These are close to, but somewhat steeper than the me- 
dian values of the {V — R), (F — /), and {R — I) gradients 
of the DS galaxies, that are respectively -0.09, -0.08, -0.03 
mag/dex per decade. A Kolmogorov-Smirnov (KS) test 
shows that only the (B-V), (R-I) pair of distribution is 
marginally different. 



In the following we will investigate quantitatively the 
evolutions of color gradients under different assumptions 
for their origin and compare the models to the observed 
distributions. 



4. The Models 
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Table 6. The color gradients of the Local Sample. Column 1 gives the galaxy name, Col. 2 the galaxy cluster (F for 
Fornax, V for Virgo, C for Coma), Col. 3 the U-B logarithmic gradient, Col. 4 its error. Col. 5 the U-B color at i?e/2. 
Col. 6 its error. Col. 7 the source. Col. 8 the U-V logarithmic gradient. Col. 9 its error. Col. 10 the U-V color at i?e/2, 
Col. 11 its error. Col. 12 the source. Col. 13 the B-V logarithmic gradient. Col. 14 its error. Col. 15 the B-V color at 
i?e/2. Col. 16 its error. Col. 17 the source. The source codes are: F for Franx et al. 1989, G for Goudfrooij et al. 1994, 



P for Peletier et al. 1990, 1 indicates that the V profile was obtained interpolating the B and R profiles (see text). 



Errors are due to sky subtraction and therefore systematic. 



Galaxy CI a(U-B) da U-B dCol S a(U-V) 



da U-V dCol 



i(B-V) 



da B-V dCol 



A 496 


- 


-0.05 


0.10 


0.57 


0.01 


p 


-0.28 


0.11 


1.45 


0.12 


I 


-0.10 


0.02 


0.94 


0.03 


I 


IC 1101 


- 


-0.05 


0.10 


0.69 


0.01 


p 


-0.09 


0.12 


1.71 


0.15 


I 


-0.06 


0.03 


0.97 


0.04 


I 


NGC 1379 


F 


-0.06 


0.18 


0.40 


0.09 


F 


-0.02 


0.28 


1.35 


0.03 


I 


-0.05 


0.04 


0.92 


0.02 


I 


NGC 1399 


F 


-0.08 


0.15 


0.55 


0.07 


F 


-0.15 


0.12 


1.52 


0.06 


F+G 


-0.07 


0.01 


0.97 


0.02 


G 


NGC 1404 


F 


-0.07 


0.04 


0.57 


0.01 


F 


-0.08 


0.03 


1.53 


0.01 


F+G 


-0.01 


0.01 


0.96 


0.01 


G 


NGC 1427 


F 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


-0.05 


0.01 


0.89 


0.01 


G 


NGC 4365 


V 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


-0.04 


0.01 


0.97 


0.01 


G 


NGC 4374 


V 


-0.17 


0.02 


0.54 


0.01 


P 


-0.22 


0.04 


1.57 


0.02 


I 


-0.04 


0.01 


1.02 


0.01 


I 


NGC 4387 


V 


-0.01 


0.02 


0.49 


0.01 


P 


0.01 


0.08 


1.54 


0.02 


I 


-0.02 


0.02 


1.05 


0.01 


I 


NGC 4406 


V 


-0.15 


0.02 


0.52 


0.01 


P 


-0.16 


0.03 


1.52 


0.02 


I 


-0.02 


0.01 


1.00 


0.01 


I 


NGC 4472 


V 


-0.15 


0.03 


0.55 


0.01 


P 


-0.18 


0.05 


1.58 


0.03 


I 


-0.03 


0.01 


1.05 


0.01 


I 


NGC 4478 


V 


-0.12 


0.04 


0.52 


0.01 


P 


-0.14 


0.07 


1.52 


0.02 


I 


-0.03 


0.02 


0.99 


0.01 


I 


NGC 4486 


V 


-0.23 


0.04 


0.41 


0.01 


P 


-0.27 


0.03 


1.37 


0.01 


P-hG 


-0.04 


0.01 


0.96 


0.01 


G 


NGC 4564 


V 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


-0.12 


0.01 


0.93 


0.01 


G 


NGC 4621 


V 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


-0.03 


0.01 


0.92 


0.01 


G 


NGC 4636 


V 


-0.22 


0.04 


0.46 


0.01 


P 


-0.27 


0.07 


1.45 


0.05 


I 


-0.05 


0.02 


0.99 


0.01 


I 


NGC 4649 


V 


-0.15 


0.04 


0.65 


0.01 


P 


-0.18 


0.06 


1.69 


0.03 


I 


-0.04 


0.01 


1.04 


0.01 


I 


NGC 4660 


V 


- 


- 


- 


- 


- 


- 


- 


- 


- 


- 


-0.05 


0.01 


0.98 


0.01 


G 


NGC 4874 


c 


-0.07 


0.06 


0.55 


0.01 


P 


-0.20 


0.05 


1.50 


0.05 


I 


-0.10 


0.02 


0.97 


0.01 


I 


NGC 4889 


c 


-0.08 


0.04 


0.72 


0.01 


P 


-0.11 


0.09 


1.65 


0.04 


I 


-0.10 


0.03 


0.95 


0.02 


I 


NGC 6086 


- 


-0.15 


0.04 


0.62 


0.01 


P 


-0.17 


0.14 


1.58 


0.06 


I 


-0.03 


0.03 


0.96 


0.01 


I 



4-.1. The Basic Tool: Simple Stellar Populations 



Our modeling is based on the assumption that the light- 
averaged stellar population in the Es between z « 0.4 
and z — can be described by Simple Stellar Populations 
(SSP), i.e. coeval and chemically homogeneous assembhes 
of single stars. The SSP models for this work are computed 
with the evolutionary synthesis code described in Maras- 
ton (1998), in which the Fuel Consumption theorem (Ren- 
zini & Buzzoni 1986) is used to evaluate the energetics of 
Post Main Sequence stars. The input stellar tracks for the 
SSP models used here are from Cassisi {private communi- 
cation; see also Bono et al. 1997), covering the metallicity 
range [Fe/H] = [—1.35-^0.35], with helium-enrichment pa- 
rameter AY/AZ=2.5. The amount of mass loss along the 
Red Giant branch (RGB) is parametrized a' la Reimers 
(1975), with the efficiency 77=0.33 as calibrated by Fusi 
Pecci & Renzini ( 1976| ). The mass loss prescriptions af- 
fect the evolutionary mass on the horizontal branch (HB), 
leading to warmer HBs as the mass loss increases. We 



adopt Salpeter IMF (*(to) = Am-'-^+''\ a;=1.35) down 
to a lower mass limit of 0.1 M©. For plausible variation of 
the IMF slope, the opti cal co lors of SSPs are virtually un- 
changed (e.g. Maraston 1998). The evolutionary synthesis 



code has been updated for the computation of the SSP 
Spectral Energy Distributions (SEDs), as a function of t 
and [Fe/H]. We adopted the spectral library of Lejeune, 
Cuisinier and Buser (1998) to describe the stellar spec- 
tra as functions of gravity, temperature and metallicity. 
A redshifted grid of SSP model SEDs has also been com- 
puted to interpret the colors of the DS galaxies. We have 
applied the same redshift z = 0.3775 to all the models. 
This is the average value of the z distribution of the DS 
sample, having assigned to the galaxies without redshift 
determination, the average redshift value of the objects 
with a measured redshift. The uncertainties introduced by 
a ~ 0.03 variation in redshift are of the same order as the 
statistical errors. Model SEDs have been convolved with 



the UBVRI Johnson-Cousins filter functions (Buser 1978 



Bessel 1979) to yield synthetic magnitudes. The applied 
zero points are (C/,B,y,i?,/)=(0.02, 0.02, 0.03, 0.0039, 



0.0035) for the Vega model atmosphere of Kurucz ( |1979| ) 
with (Te, g, [Fe/H])=(9400, 3.9, -0.5). This ensures consis- 
tency with our HST photometry (Holtzman et al. 1995). 
We refer to Maraston (EOOO) for more details. 

Figures 1^ and ITol summarize the properties of the mod- 
els relevant to our case. Fig. g shows the {U — B), {U — V) 
and (B — V) colors of the rest-frame SSP models, together 
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Fig. 6. Top-left: the color gradients of the distant sam- 
ple as a function of i?e(675) in arcsec. Open squares refer 
to the {V — R), crosses to the {V — I) and filled squares 
to the {R — /) colors. The bars show the statistical er- 
rors. Bottom-left: the color gradients of the distant sample 
as a function of tot (675). The bars show the systematic 
errors. Top-middle: the cumulative distribution of ratio 
between the gradient and the statistical error. The solid 
line refers to the {V — R), the dotted to the {V — I), the 
dashed to the (R — I) colors. The vertical line shows the 3cr 
limit. Bottom-middle: the correlation between statistical 
and systematic errors. Symbols as above. Top-right: the 
correlation between the {V~R) and the {R~I) gradients. 
Statistical errors are shown. Bottom- right: the correlation 
between the {V — T) and the {R — I) gradients. Statistical 
errors are shown. 



with their derivatives with respect to age and metallicity 
(given as [Fe/H]). Derivatives are computed as finite dif- 
ferences on the grid and assigned to the midpoints in age 
and metallicity. 

In general, optical colors become redder with increas- 
ing age and/or increasing metallicity. In our models, there 
are two exceptions to this trend, both at low metallicities. 
The [Fe/H] = -1.35 models become bluer in U — B as they 
age from ~ 3 to '^ 5 Gyr. This is due to the properties 
of thcU -B stellar colors for T^ - 8500-7500 K (see e.g. 



Castelli 1999 ) typical of turnoff stars in this age range. At 
higher metallicities this effect does not appear, because 
when the turnoff stars have temperatures in this range, 
their contribution to the total light is less important due 
to the presence of AGB stars (Maraston [2000 ). 

At old ages, the same set of models become bluer and 
bluer in B — V from ~ 12 Gyr onward, reflecting warmer 



Fig. 7. Cumulative distributions of the B-band absolute 
magnitudes of the z=0 (solid) and z sa 0.4 (dotted) sam- 
ples. 



.. V-R(az = 0.37f 

: U-B@z = 

: KS = 0.567 



.. V-IiS>z = 0-37, 
.iT _ U-V@z = 
:; KS = 0.248/ 



.. R-I@z = 0.37i 
.•■Jl _ B-V@z = ) 

:: KS - 0.037 ■ 



-0.3 -0.2 -0 



-0.3 -0.2 -0.1 
Gradients 



-0.1 -0.05 



Fig. 8. Cumulative distributions of color gradients of the 
z=0 (solid) and z « 0.4 (dotted) samples. The left, central 
and right panels show respectively the {U — B), {U — V) 
and {B — V) gradients of the LS, to be compared to the 
{V -R), (y - I) and {R - I) gradients of the DS. Each 
panel is labeled with the results of the KS tests applied to 
the corresponding pairs of distributions. 



and warmer HBs. With r; = 0.33, this effect sets in for 
[Fe/H]^ —0.5. It disappears when no mass loss is ap- 

The 



1988) 



plied (Maraston 2000 see also Chiosi et al. 
time derivative of colors shows a noisy behavior due to 
the discrete modeling. However, a general trend can be 
appreciated, with the derivative initially decreasing and 
leveling off at ~ 0.01 - 0.04 mag/Gyr for ages older than 
~ 7 Gyr. The time derivatives appear almost independent 
on metallicity for Z > 0.5 Zq, while the lowest Z set of 
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models behave differently due to the two effects described 
above. 

The derivatives of colors with respect to metallicity 
show a different trend with age, with colors generally be- 
coming more and more sensitive to [Fe/H] as the SSPs 
get older. It is worth noticing that the plotted curves are 
color variations per metallicity decade. The dependence 
on Z is quite different, with low metallicity model col- 
ors much more sensitive to Z variations than the high Z 
ones. We preferred to plot the derivatives with respect to 
[Fe/H] because these are the actual tools used in our sim- 
ulations. We also point out that in our models the helium 
abundance varies in lock-steps with Z. This influences the 
temporal behavior of the SSP colors and of their deriva- 
tives, especially at high metallicity. 

Figure |ll| compares our models (shown down to 
you nger ages for the purpose) to th e ones of Vazde kis et 
al. ( |1996[ ), Bruzual & Chariot ( |l996[ ) and Worthey ( [1994D 



11998 
1996| at 



as a function of metallicity and age (see Maraston 
for a comparison to the models of Tantalo et al 
solar metallicity). On the whole, the agreement between 
all the models in the {B — V^ color is good. The variations 
between the models are of the same order of the uncertain- 
ties in the observed colors. Our SSP optical colors compare 
remarkably well with Vazdekis et al. ( 19961 ). Larger differ- 
ences between the models are observed in the {U — B). 



Bruzual & Chariot (1996) models tend to have system- 
atically redder {U — B). Worthey (1994) models produce 
instead systematically bluer (U — B) (and [U — V)) colors 
at a given {B — V). A thorough investigation of these dis- 
crepancies will be presented in Maraston (pOOO). The age 
(at fixed metallicity) or metallicity (at fixed age) differ- 
ences inferred using the different models between the cen- 
tral and outer parts of the galaxies (the galaxy NGC 1399 
is shown as an example) are consistent within the errors. 
The results on the origin of the color gradients presented 
here also depend on the derivatives of the colors with re- 
spect to age and metallicity. This dependence is difficult 
to estimate quantitatively because of the coarseness of the 
grid of the other sets of models, and of the unavailability 
of their entire SEDs. The close similarity of our models to 
the Vazdekis is reassuring. 

Similar to Fig. 0, Fig. O shows the properties of the 
redshifted SSP SEDs in the V, R and / bands. Broadly 
speaking, the {V — R), {V — I) and (i? — /) colors trends 
with age resemble those of the {U—B), {U—V) and {B—V) 
colors. Some differences are present, as a result of the dif- 
ferent effective wavelengths and sensitivity curves of the 
corresponding pairs of filters. In particular we notice that 
the time derivative of the {V — I) color is approximately a 
factor of 2 smaller than that of the [U — V) color, reminis- 
cent of Fig. g. Thus we expect that the (JJ — V) colors of 
nearby galaxies are more sensitive to age gradients than 
the {V — I) colors of distant objects. 

Simple inspection of Figs. |^ and |l^ already allows us 
to define the framework for the interpretation of color gra- 



dients in ellipticals. The filled and open dots show the 0.2 
and 2 Rf, colors of the LS and DS galaxies. It is clear that 
the range of ages and metallicities considered here is ad- 
equate to model the observed colors both as a pure age 
sequence (at a fixed metallicity) or as a pure metallicity 
sequence (at a fixed age) . 

Turning now to the gradients' evolution, we start con- 
sidering that the LS galaxies all (with one exception) have 
negative gradients, with colors being redder in the central 
regions. Thus, if they are caused by age gradients only, 
centers must be older than the outer parts (~ 5 Gyr per 
radial decade). In addition, if the present age of the Es 
centers is large enough (« 15 Gyr), we expect the age 
derivatives of the colors to stay approximately constant 
for look-back times up to « 7 Gyr, and therefore to mea- 
sure color gradients in the DS of approximately the same 
size in (V — R) and {R — I), and twice as small in {V — I). 

In contrast, if color gradients are caused by metallicity 
gradients only, the centers of local Es must be more metal 
rich than the outer parts (« 1.6 more metal rich per radial 
decade). In addition, given the increase with age of the 
metallicity derivatives of the colors, we expect the color 
gradients of the DS to be fiatter than those of the LS. This 
effect should be less prominent if the centers of ellipticals 
are old. 

As noted in Sect. 0, some of the gradients measured in 
CL0949 Es are positive. In Sect. 6.3 we shall argue that 
these are caused by measurement errors, with an under- 
lying distribution of intrinsically negative gradients. Here 
we note that a change in sign with time of the color time 
derivative of SSP models is possible only at low metallici- 
ties and/or low ages. Color gradients caused by metallicity 
cannot change sign with time following passive evolution. 

4-. 2. Modeling the Evolution of Color Gradients 

The detailed modeling of the color gradients evolution 
of each galaxy involves a number of steps. As a start- 
ing point, a cosmological model should be chosen to fix 
the look-back time at z « 0.4. Here we adopt Ho=65 
Km • s-i • Mpc-i and 17^=0.2, with A = 0.8. The age 
of the Universe is w 16.18 Gyr and <ib = 4.53 Gyr at 
2=0. 3775. The choice of the cosmological parameters is 
not crucial to model color gradients, but a much younger 
universe (for example, 10 Gyr) would strongly disfavor 
the age hypothesis for the origin of color gradients (see 
below). In addition, the look-back time has an influence 
on the exact value of the predicted colors. As we will show 
later, the optimal agreement between the observed color 
distributions in the LS and the DS samples is achieved for 
larger look-back times, i.e., iib ~ 8 Gyr. 

Inspection of Figs. pj-H shows that the color gradients of 
the DS galaxies are measured in the radial range 0.2—2 Rg. 
The same applies to most of the local galaxies. Therefore, 
we choose this radial decade to model the gradients. A 
more innerly defined decade (i.e., 0.1 — 1 -Re) suffers of in- 
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Fig. 9. The colors, time and [Fe/H] derivatives of the rest-frame SSP models in the U, B and V bands. The different 
line type of the color and time derivative rows refer to the metallicitics indicated in the top- left plot. The line types 
and their corresponding metallicities (midpoints of the above) of the [Fe/H] derivative are given in the bottom panel of 
the central figure. The filled dots show the 0.2 Rg colors of the LS galaxies with ages assigned on the solar metallicity 
line. The open dots show their 2 R^ colors at the (arbitrary) age of 14 Gyr. 



sufficient resolution for the more distant sample. A more 
outerly defined decade (i.e., 0.3 — 3 Re) is poorly investi- 
gated in the local sample. Thus, we analyse a given color 
gradient by examining the colors at the 0.2 Re (hereafter, 
the central region) and 2 Re (hereafter, the outer region) 



isophotes, as derived from the values listed in the Table 



We need now to choose representative SSP models for 
the inner and outer regions. Due to the age-metallicity 
degeneracy, one parameter has to be fixed to derive the 
other from the observed colors. For the central stellar pop- 
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Fig. 10. The colors, time and [Fe/H] derivatives of the redshiftcd SSP models in the V, R, I bands. The different Une 
type of the color and time derivative rows refer to the mctallicitics indicated in the top-middle plot. The line types 
and their corresponding metallicities (midpoints of the above) of the [Fe/H] derivative are given in the bottom-left 
plot. The vertical hnes show the age (5.47 and 10.47 Gyr) of galaxies at z=0.3775 that today are 10 and 15 Gyr old, 
respectively. The filled dots show the 0.2 Re colors of DS galaxies with ages assigned on the solar metallicity line. The 
open dots show their 2 Re colors at the (arbitrary) age of 14 Gyr. 



ulations we consider two options: 1) Gase m, in which all 
the galaxy centers have the same age and the colors are 
used to infer their metallicities; 2) Case a, in which all the 
galaxy centers have the same metallicity, and their colors 
are used to infer their ages. In Case m we explore two 



possible current ages: 10 and 15 Gyr; in Gase a we as- 
sign solar metallicity to all the galaxy centers. A smaller 
metallicity (i.e., half-solar) would imply implausibly large 
(larger than the age of the Universe) ages for the galax- 
ies (see Figs. U and [l0|). The canonical interpretation of 
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Fig. 11. The (U — B) (t op) an d {B — V) rest-frame color time evolution of our model (full lines), Vazdekis et al. ( 1996 , 
dotted lines), Worthey ( 1994 , short-dashed lines), Bruzual & Chariot ( 1996 , long-dashed lines) as a function of age 
(for three metallicities) and metallicity (15 Gyr age for our models and Bruzual & Chariot, 17 Gyr for Vazdekis and 
Worthey's models). The error bars show the allowed range of the inner (at 0.2i?e) colors of the LS galaxies derived 
from the gradient fits. The filled circles and squares show the inner and outer {2Re) colors, respectively, of the galaxy 
NGC 1399. 



the Color-Magnitude relation of cluster early type galax- 
ies is consistent with Case m (e.g. Arimoto & Yoshii 1986, 
Kodama & Arimoto 1997); Case a is nevertheless worth 
exploring. 



As for the gradients, we consider the two extreme pos- 
sibilities that they are driven by pure age (A) or pure 
metallicity (M) gradients. Thus we have four scenarios in 
which we model the evolution of the color gradients: 
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aA: age drives the color gradients. Central and outer SSPs 
have solar Z and colors reflect their ages. The gradi- 
ents are expected to steepen with look-back time if the 
centers are young; to remain nearly constant if they are 
old enough. 

mA: age drives the color gradients. The central SSPs of local 
Es are 10 or 15 Gyr old (mAlO and mA15); the outer 
SSPs have the same metallicity as the central ones, but 
different ages. The gradients' evolution is expected to 
be similar to case aA. 

aM: metallicity drives the color gradients. Central SSPs 
have solar Z\ the outer SSPs have the same age as 
the centers, but different metallicities. Color gradients 
are expected to flatten with look-back time. 
mM: metallicity drives the color gradients. The central and 
outer SSPs of local Es are coeval and 10 or 15 Gyr 
old (mMlO and mM15), and the colors reflect their 
metallicities. The gradients' evolution is expected to 
be similar to case aM. 

To summarize, we model the evolution of the color 
gradients of each galaxy in the LS in the four scenarios 
by: 

1.) fixing the age and metallicity of the central and outer 

SSPs according to the specific option (aA,aM,mA and 

mM); 
2.) rejuvenating the SSPs by the adopted look-back time 

(4.53 Gyr); 
3.) computing the color gradients of the rejuvenated and 

redshiftcd SSPs in the V, R and I bands. 

For consistency, we also perform the reverse experi- 
ment and compute the predicted U,B,V color gradients 
of the CL0949 galaxies after t^B (i-e., 4.53 Gyr) of pas- 
sive evolution. We tested that the results we are going 
to discuss remain the same within the errors if the whole 
color profiles are modeled at each radial distance, and a 
logarithmic gradient is fit to the predicted profiles at z=0. 

5. Ages and Metallicities of the Galaxies in the 
Four Scenarios 

We consider now the applicability of the four scenarios to 
our local and distant galaxy samples. Indeed, the ages or 
metallicities derived from the three colors available should 
agree within the measurement errors. The distribution of 
ages derived for cases aA, aM and mA for the LS and the 
DS galaxies should be similar, given the look back time. 
The distributions of metallicities derived for cases aM, 
mA and mM should also be similar. 



samples in both regions. The ages derived from different 
colors are in agreement (within the rather large statistical 
errors), better for the DS objects. We note that the LS 
objects appear in the plot only when the specific pair of 
colors is available. However, the average ages and rms in- 
dicated in the plots are those obtained considering all the 
objects with the specified color available. On the average, 
the centers of LS galaxies appear ~ 12 Gyr old, while the 
peripheries are ^ 7 Gyr old. The DS galaxies appear ^ 4 
Gyr old in the center, ~ 2.5 Gyr old in the peripheries. 
Thus, the consistency of the aA hypothesis is rather ques- 
tionable, both because the average age difference between 
the centers of the two samples is larger than our adopted 
look-back time, and because the age gradients implied un- 
der this scenario are much smaller for the DS sample than 
for the LS one. In addition, the colors of some local galax- 
ies require a present age younger than 4.5 Gyr or an age 
older than the age of the universe in our chosen cosmol- 
ogy. In other words, only a subsample of the LS galaxies 
could be the passively evolved descendants of the CL0949 
objects. This is partly a consequence of forcing a solar Z 
for all the considered SSPs. 




U-B 



Fig. 12. Ages assigned to the center (0.2 Re) and the outer 
isophotes (2 Re) of the local galaxies, assuming that they 
have solar metallicity. The given mean age and rms are 
derived using all the galaxies available in a given filter. 
The cross shows the mean of the measurement errors. 



5.1. Age gradients: the aA Case 



Figures 12 and [13 show the ages assigned to the center 



(0.2 Re) and the outer isophotes (2 Re) of the local and 
distant galaxies, having assumed solar metallicity for both 



5.2. Metallicity gradients: the mM. Case 

Figures |lj and |l5| show the metallicities assigned to the 
center (0.2 Re) and the outer isophotes (2 Re) of the local 
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Fig. 13. Ages assigned to the center (0.2 Re) and the outer 
isophotcs (2 Re) of the distant galaxies, assuming that 
they have solar metallicity. Labels and crosses as in Fig. 

HI 



and distant galaxies, assuming a present day age of 15 Gyr 
for the LS objects. At z = the metallicities derived from 
the three colors agree well. Local galaxies appear to be 
slightly subsolar at their centers if 15 Gyr old, solar if 10 
Gyr old. The central SSPs of the DS objects require only 
half solar metallicity, if the look-back time is 4.53 Gyr, 
but with tjb ~ 8 Gyr the central Z is larger by 0.15 dex. 
This discrepancy is the same present in the aA case: for 
our adopted look-back time the observed central colors of 
the DS objects are too blue (those of the LS are too red) 
when compared to passively evolving models. The metal- 
licity gradients are « 0.2 dex for both DS and LS, indi- 
cating that metallicity driven color gradients are a viable 
interpretation (see Sect. H). However, we notice that the 
DS galaxies span a wider range in metallicity gradients, 
compared to LS objects. This indicates a slight inconsis- 
tency also for the mM scenario, but see Sect. 6.3 for a 
possible way out. 



5.3. The aM. and mA Cases 

The aM case is always applicable: since the central ages 
derived in ^^ are ~ 12 Gyr on the average, and the central 
metallicities are supposed to be solar, the other colors for 
both DS and LS galaxies can be reproduced consistently 
in all filters by decreasing the metallicity by 0.2 dex. 

The applicability of the mA scenario strongly depends 
on the assumed central age. When 10 Gyr is adopted for 
the centers of LS galaxies, a substantial number of objects 



Fig. 14. Metallicities assigned to the center (0.2 Re) and 
the outer isophotes (2 Re) of the local galaxies, assum- 
ing that their age is 15 Gyr. The given mean metallicity 
and rms are derived using all the galaxies available in a 
given filter. The cross shows the mean of the measurement 
errors. 
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Fig. 15. Metallicities assigned to the center (0.2 Re) and 
the outer isophotes (2 Re) of the distant galaxies, assum- 
ing that their age is 10.47 Gyr. Labels and crosses as in 
Fig. H. 
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appear to have outer SSPs younger than our adopted look- 
back time. In addition, the average age of the outer regions 
of the DS galaxies is « 3 Gyr, too old to evolve into the 
LS galaxies. The situation is different if the central age of 
local galaxies is fixed to 15 Gyr. The average outer age of 
the LS is « 10 Gyr in this case, with only a few objects 
younger than our adopted look-back time. The peripheries 
of the DS galaxies appear w 6 Gyr old on the average, 
which is consistent with the average age of the LS. 

To summarize, the metallicity driven color gradients 
hypothesis is always applicable, while the age hypothesis 
holds only if local galaxies are old (~ 15 Gyr). This re- 
flects the fact that a small metallicity variation is required 
to explain the observed color gradients, while a large age 
gradient (of ^ 5 Gyr) is needed to fit the color gradients of 
the LS objects. With a 4.5 Gyr look-back time, a substan- 
tial number of LS galaxies do not have their corresponding 
object in the DS. 

6. Results 

In this section we present the results of the detailed mod- 
eling of the evolution of the central colors and color gra- 
dients under the four scenarios, according to the proce- 
dure outlined in Sect 



4.2 



We will compare the model 
predictions to the data by considering the cumulative dis- 
tribution functions of the color gradients and applying 
Kolmogorov-Smirnov tests to the pairs of correspondent 
theoretical and observational distributions. In particular, 
the redshifted V, R, I colors of rejuvenated LS galaxies 
will be compared to those of CL0949 objects; the U, B, V 
colors of passively evolved DS galaxies will be compared 
to the U, B, V colors of local ellipticals. As discussed 
above, the colors of some LS galaxies are so blue that 
ages shorter than Ilb would be required. In these cases 
the object "drops out" of the sample by getting an as- 
signed (blue) color of -1. For some choices of age or metal- 
licity the colors are outside our grid of models. In these 
cases we assign the minimum (-1.35) or maximum (0.35) 
metallicities, and 30 Myr of age (the smallest age of the 
grid), if the extrapolated age is smaller than 0, or 15 Gyr, 
if larger ages are extrapolated. 

6.1. The distributions of central colors 

Fig. nq compares the cumulative distributions of observed 
and simulated central colors of local and distant galax- 
ies. The Kolmogorov-Smirnov probability KS that the 
observed and simulated data are drawn from the same 
underlying distribution are always smaller or equal than 
0.001. None of the simulated central color distributions 
appear compatible with the observed ones. This result is 
confirmed using the Wilcoxon test. The simulated cen- 
ters of the DS objects appear too blue (by « 0.04 mag in 
B — V) when compared to the local ones, and, vice versa, 
the centers of local galaxies appear too red (~ 0.03 mag 



in R-I) when compared to the CL0949 galaxies. The dif- 
ferences are more pronounced in the {U — B) and {V — R) 
and {U — V) and {V — /) pairs. Measurement errors alone 
seem insufficient to explain the difference between these 
pairs (see |6.3D , while the (B — V) and {R — I) colors might 
not differ significantly in a statistically sense. The uncer- 
tainties in the galaxy redshifts and in the calibration to 



the Johnson-Cousins system (see Holtzmann et al. 1995) 



combined with the notoriously difficult photometric cal- 
ibration of the U band probably suffice to explain our 
findings. In any case, we tested that the conclusions con- 
cerning the origin of color gradients are not affected by 
this possible inconsistency of the central colors. The dis- 
tributions of gradients shown in Fig. |l^ are hardly changed 
when the look-back time is increased from 4.53 to 8 Gyr. 
Assuming that the color differences are real, they can 
be explained using a larger look-back time. For example, 
tLB = 8 Gyr (obtained with the same cosmology and a 
rather extreme low value of the Hubble constant, Hq ~ 37) 
gives reasonable KS probabilities (5-80% for the different 
colors and models). The result might also point to recent, 
minor episodes of star formation in the distant galaxies 
(Poggianti et al. 1998, Couch et al. 
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Fig. 16. The cumulative distributions of observed (solid 
line) and simulated (for t^s — 4.53 Gyr) central colors of 
z=0 (top panels) and z « 0.4 (bottom panels) galaxies. 
The dotted lines show Case a, the short- and long-dashed 
line the Case m for present-day central ages of 10 and 
15 Gyr. The Kolmogorov-Smirnov probabilities that the 
observed and simulated data are drawn from the same 
underlying distribution arc always smaller or equal than 
0.001. 
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6. 2. The distributions of color gradients 

Fig. ^ compares the cumulative distributions of the ob- 
served and simulated color gradients. The Kolmogorov- 
Smirnov probabilities that the observed and simulated 
data are drawn from the same underlying distribution arc 
given in Table 0- 

The a A Case is clearly ruled out. The very small dif- 
ferences between the central and outer ages required at 
z w 0.4 (see Fig. 03) evolve into too shallow gradients 
at z=0. Along the same line, the 5 Gyr age difference in- 
ferred at z=0 (see Fig. n2) produces too steep gradients 
at z « 0.4, because the color time derivatives incr ease for 
decreasing age (see Fig. |y and discussion in Sect. 4.1). 



The mM case is viable for both the assumed present- 
day central ages of 10 and 15 Gyr. The t-test confirms 
that the means of the observed and simulated distribu- 
tions in all colors at both z==0 and z w 0.4 are statistically 
identical. The f-test on the rms of the distributions con- 
firms the slightly lower KS probabilities obtained for the 
B — V gradients for z=0 and for the (B. — I) for z« 0.4. 
As discussed in Sect. 4T, all (except one) local gradients 
are negative, while some positive gradients are observed 
at z « 0.4. This makes the observed widths of the z « 0.4 
gradient distributions larger than the simulated ones. Er- 
rors might partly explain this effect (sec discussion in Sect. 
6.3). Indeed, the widths of the distributions are somewhat 
dissimilar, with the DS presenting a tail of steep gradients 
with no counterpart in the simulated LS. Since the metal- 
licity color derivatives decrease with decreasing age, the 
evolution leads to shallower (not steeper) gradients. The 
same considerations apply to the aM Case, which provides 
also an acceptable model of color gradients. 

The mA Case does not offer a perfect clear-cut solu- 
tion. One can clearly exclude the case of central present- 
day young ages (10 Gyr), since the resulting outer ages at 
z=0 are either younger than the assumed t^B , or so young 
to correspond to a very steep rejuvenated color gradient. 
The outer colors of the distant galaxies require such large 
ages that too shallow gradients are produced at z=0. In 
contrast, present day central ages of 15 Gyr still make the 
age hypothesis a possible explanation for color gradients. 
In this case, the outer ages (some « 5 Gyr smaller than 
the central ones) required to match the colors of the z=0 
galaxies are large enough (see Sect. 5.3) to produce mild 
gradients at z « 0.4, even if a tail of too steep gradients 
is produced, which is not observed in the data. Similarly, 
only a small fraction (w 10 — 20%) of the distant galaxies 
cannot be matched in their outer colors by this mecha- 
nism. The remaining objects are assigned outer ages that 
produce reasonable gradients locally. 

6.3. The influence of errors 

In the previous sections we pointed out that some of the 
problems recognized in the modeling (the difference in the 
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Fig. 17. The cumulative distributions of observed and 
simulated color gradients of z=0 (top panels) and z « 0.4 
(bottom panels) galaxies. In each of the plots the thick 
solid line shows the observed distributions. Two rows of 
plots are shown for the local and distant case. The top row 
refers to the age-driven color gradients. Here the dotted 
lines show the aA case, the short-dashed lines the mAlO 
case, the long-dashed lines the mA15 case. The bottom 
row refers to the metallicity-driven color gradients. Here 
the dotted lines show the aM case, the short-dashed lines 
the mMlO case, the long-dashed lines the mM15 case. 
The Kolmogorov-Smirnov probabilities that the observed 
and simulated data are drawn from the same underlying 
distribution are given in Table M. 

median central colors of the LS and DS samples, the pres- 
ence of positive gradients at z « 0.4, the differences in 
the observed and predicted width of the gradient distribu- 
tions) might be explained taking into account the observa- 
tional errors. Errors are dominated by systematics for the 
LS galaxies, and have both a statistical and a systematic 
component for the DS. They affect both the observed and 
the predicted distributions. We estimate the effect of sta- 
tistical errors on the observed distributions by bootstrap- 
ping the measured colors or color gradients with their sta- 
tistical errors, assumed to be gaussianly distributed. The 
estimation of the effect on the predicted color and color 
gradient distributions requires the choice of a model. For 
simplicity, we focus on the mM Case with present-day age 
of 10 Gyr. We bootstrap the measured colors at 0.2 and 
2 Be with their statistical errors at z « 0.4, assumed to 
be gaussianly distributed. We compute the metallicities 
implied by the SSP models, predict the expected colors at 
z=0, and derive the color gradients and their cumulative 
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Table 7. The Kolmogorov-Smirnov probabilities of the color gradient distributions. Col. 1 gives the examined Case, 
Cols.2-5 the KS probabihties for the {U - B), (U - V), (B - V) color gradients at z = 0, Cols. 6-7 for the {V - R), 
(V -I), (R-I) at z«0.4. 



Case 


KSu-B 


KSu-v 


KSb-v 


KSv-R 


KSv^i 


KSr-i 


aA 


7.5e-05 


4.4e-06 


2.9e-07 


0.05 


0.013 


0.002 


mMlO 


0.79 


0.50 


0.18 


0.62 


0.50 


0.18 


mM15 


0.67 


0.51 


0.09 


0.17 


0.29 


0.22 


aM 


0.37 


0.15 


0.20 


0.83 


0.50 


0.32 


mAlO 


0.016 


0.004 


0.08 


0.009 


0.0002 


2.7e-05 


mA15 


0.83 


0.50 


0.01 


0.62 


0.87 


0.18 



distributions. Figs. [T| and |T^ show the result of 30 such 
bootstraps. The typical width of the obtained measured 
and predicted color distribution clouds is « 0.02 mag. As 
discussed in Sect. 0, the statistical errors dominate over 
the systematic effects. The typical width of the obtained 
measured and predicted distributions clouds of the color 
gradients is « 0.02 — 0.04 mag/dex. In this case the sys- 
tematic effects bracket the statistical variations. 
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Fig. 18. The influence of statistical and systematic errors 
on the cumulative distributions of the observed and pre- 
dicted central colors for the LS (upper) and the DS (lower) 
galaxies. For the DS galaxies and for each color we show 30 
distributions (dotted lines) obtained by bootstrapping the 
measured colors and modeling their evolution (see text for 
explanation) . The thick lines show the measured distribu- 
tions. The dashed thick lines show the effects of system- 
atic errors. For the LS galaxies we show the measured and 
evolved distributions taking into account the systematic 
errors (lines with crosses). 
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Fig. 19. The influence of statistical and systematic er- 
rors on the cumulative distributions of the observed and 
predicted color gradients for the LS (upper) and the DS 
(lower) galaxies. For the DS galaxies and for each color 
we show 30 distributions (dotted lines) obtained by boot- 
strapping the measured colors and modeling their evolu- 
tion (see text for explanation). The thick lines show the 
measured distributions. The dashed thick lines show the 
effects of systematic errors. For the LS galaxies we show 
the measured and evolved distributions taking into ac- 
count the systematic errors (lines with crosses). 



Fig. n8^ shows that measurement errors alone are insuf- 
ficient to explain the difference in central color between the 
LS and DS samples, because the simulated and observed 
distributions at both z=0 and z « 0.4 are separated even 
when errors on both sides are considered. A larger look- 
back time (8 Gyr) produces a good overlap. 

The inspection of Figure [l9| shows first that the re- 
sults described in Sect. 3.2 are robust against the errors. 
Cases a A and mA with a central present age of 10 Gyr 
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are ruled out with high significance because the large tail 
of steep negative gradients that they produce at 2; w 0.4 
is not present in any of the distributions with errors. As 
a second remark, errors help explaining the discrepancies 
between the observed and modeled distributions. Com- 
plete overlap at both z=0 and z « 0.4 is achieved between 
the observed and modeled distributions of the U-B and 
V — R gradients. The agreement for the U — V and V — I 
gradient distributions is improved, with most of the differ- 
ences in the widths accounted by the errors. The 10% tail 
of distant galaxies with the largest positive R-I gradients 
(objects 211 and 312) cannot be explained by the observa- 
tional errors. However, the two galaxies are the smallest 
of the sample and it can be argued that their gradients 
are not well determined (see Sect. 0). Therefore, we con- 
clude that the mM, aM, and, to less extent, mA15 cases 
provide a satisfactory description of the color gradients of 
(bright) cluster ellipticals at z=0 and z « 0.4 following 
passive evolution. 

6.4- The role of dust 



Wise & Silva (1996) examined the possibility that color 
gradients in elliptical galaxies are caused entirely by dust 
effects. They fit the color profiles of the LS galaxies con- 
sidered here adopting different radial distributions and to- 
tal mass of the dust component. Their best fitting mod- 



els are characterized by a distributions oc 



and dust 



masses much higher than those directly determined from 
the IRAS data. The ratios of these masses are in the range 
3.6-307, only NGC 1399 having ratio of order unity. There- 
fore we are confident that the radial variation of the stellar 
populations plays a major role in driving the color gradi- 
ents. Nevertheless, part of the gradient may be due to 
dust. We have thus tested the hypothesis that 50% of the 
color gradients observed are caused by dust, to verify how 
robust our conclusions of Sect. 1x2 are. We reduce the gra- 
dients by a factor 2 by making the 0.2 Re colors bluer and 
leaving the 2 R^ as measured. As a consequence, for Case 
a we require central mean ages smaller by w 2 Gyr at z=0 
but only shorter by « 0.4 Gyr at 2; ~ 0.4. For Case m we 
require central mean metallicities smaller by « 0.1 dcx. 
This does not modify much our conclusions about the vi- 
ability of the different mechanisms examined in Sect. 6.2. 
Case aA remains ruled out with high confidence. A metal- 
licity change (Cases mM and aM) of « 0.1 dex per radial 
decade is enough to match both local and distant gradient 
distributions. The Case mA with 15 Gyr central ages and 
« 2.5 Gyr younger peripheries allows a consistent model 
of all local and distant objects. The Case mA with 10 Gyr 
central ages is less unprobable than before, but still not 
acceptable. 

These findings can be used to estimate the viability 
of mixed scenarios, where combinations of age and metal- 
licity variations are driving the color gradients. The pos- 
sibility that metallicity increases in the outer regions of 



ellipticals is very unlikely, because the age driven model 
(Case aA) cannot explain color gradients already when 
the metallicity is kept constant, and works only with large 
central ages (Case mA15). In this case 0.1 dex increase in 
metallicity in the outer parts would require ~ 2.5 Gyr 
younger ages, enough to attribute ages of the order of 
tLB to 50% of the LS galaxies. The possibility that both 
metallicity and age decrease outwards can be constrained 
if the present day central ages of cluster ellipticals are < 10 
Gyr. In this case at most half of the observed color gradi- 
ents can be explained by age variations. Scenarios where 
metallicity "overdrives" the color gradients by compensat- 
ing age increases in the outer regions of the galaxies are 
also possible. Approximately 0.1 dex additional decrease 
in metallicity is needed to compensate 2.5 Gyr increase in 
age every radial decade. 

7. Conclusions 

We presented the surface brightness profiles of the 20 
brightest early-type galaxies of CL0949 at redshift z=0.35- 
0.38 from HST WF2 frames taken with the filters F555W, 
F675W, F814W. We determined the color profiles {V - 
^)('')j (^ ~ -^)(''); ^"^d (i? — /)(r), and fit logarithmic gra- 
dients in the range —0.2 to 0.1 mag per decade. These 
values are similar to what is found locally for the colors 
{U — B), {U — V), (B — V) which approximately match the 
(V-R), (V-I), (R-I) at redshift « 0.4. We analyzed the 
results with up to date stellar population models, explor- 
ing whether the following mechanisms are able to explain 
the colors of cluster ellipticals and their gradients assum- 
ing passive evolution between redshift z=0 and 0.4. Case 
aA: The central colors and their gradients are due to age 
differences. Galaxies have all solar metallicity. Case mM: 
The central colors and their gradients are due to metal- 
licity differences. Galaxies have all the same age (10 or 15 
Gyr at present). Case aM: The color gradients are due 
to metallicity. Galaxies have solar metallicity in their cen- 
ters, the central colors are defined by the age of the galaxy. 
Case mA: The color gradients are due to age. Galaxies 
have 10-15 Gyr old centers at present, their central colors 
are defined by the metallicity of the galaxy. We reached 
the following conclusions: 

— The age driven gradient with fixed-metallicity (Case 
aA) can be ruled out. It fails to reproduce the distribu- 
tion of color gradients of distant (predicting too steep 
negative gradients) and local (predicting too shallow 
gradients) galaxies. 

— The metallicity driven gradient with fixed-age (Case 
mM) is viable for both galaxy ages assumed (10 or 15 
Gyr at the present day). A metallicity gradient of « 0.2 
dex per radial decade is needed to explain the observed 
color gradients of local and distant galaxies. The gra- 
dient distributions at z=0 and 0.4 can be matched well 
by taking into account the observational errors. 
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— The metallicity driven gradients model assigning a 
fixed (solar) metallicity to the centers (Case aM) 
works also well, being very similar to case mM for 
a 10 Gyr age. 

— The age driven gradient model assigning a fixed age 
to the centers (Case mA) does not work if the present 
day central age is smaller than or equal to 10 Gyr. In 
this case too steep negative gradients at z « 0.4 are 
produced and furthermore the modeling of some z — 
galaxies is not possible, because it would require ages 
smaller than tLB- A central age of 15 Gyr allows the 
modeling of most galaxies with an age gradient of ~ 4 
Gyr per radial decade. 

— The central colors of local (bright) cluster galaxies im- 
ply central ages of « 10 Gyr (if solar metallicity is 
assumed) or slightly subsolar ([Fe/H]w —0.1) metal- 
licity, if a central age of 15 Gyr is assumed. 

— With tLB = 4.53 local galaxies appear too red at 
z « 0.4 and distant galaxies appear too blue at z=0. 
Possible reasons are the uncertainties in the galaxy 
redshifts, in the calibration to the Johnson-Cousins 
system at z w 0.4 and the U band at z = 0. The 
distributions of the central colors of local and distant 
ellipticals can also be matched by passive evolution of 
simple stellar populations if large look-back times (« 8 
Gyr) are assumed. 

— The results described are obtained consistently in the 
three pairs of colors {U ~ B and V — R, U — V and 
V -I, B-V andR-I). 

— The conclusions reached above are still valid if the 
measured gradients are reduced by a factor 2 to simu- 
late the possible presence of a screen of dust. A 0.1 dex 
reduction of metallicity per radial decade would suffice 
in Cases mM and aM to explain the data. Case mA 
with 15 Gyr present day central age and ~ 2.5 Gyr 
younger outer regions would also be a possible expla- 
nation. 

Bright cluster ellipticals appear 10-15 Gyr old at z—Q 
whatever the modeling assumptions. This is an important 
constraint that models of galaxy formation must take into 
account. 

The present results do not allow per se to discrimi- 
nate whether metallicity or age are generating the color- 
magnitude relation of cluster ellipticals. Both mM and 
aM cases provide a reasonable model of the color gradients 
observed locally and at z « 0.4 down to the magnitude 
limit of our sample {B « —21). However, the metallicity 
scenario is favored when the CM relation down to fainter 



magnitudes is considered (Kodama and Arimoto 1997). 

The presence of a metallicity gradient in galaxies (mM 
or aM case) is our favored mechanism for the origin of 
color gradients, since a fraction of galaxies cannot be mod- 



Field North and in clusters. This indicates that our result 
does not depend on our particular choice of models. 

The mild center to periphery metallicity difference (- 
0.2 dex per radial decade) reiterates the failure of classical 
monolithic collapse models with Salpeter IMF, which pre- 
dict much stronge r gra dients (~ —0.5 dex between center 
and Re, Carlberg 1984 ). In hierarchical galaxy formation 
scenarios, originally present gradients are diluted by the 
merging process (White 1980| ). This goes in the direction 
of accounting for the observations, but quantitative pre- 
dictions are still lacking. Passive evolution seems adequate 
to explain the color gradients at both z=0 and z « 0.4, 
adding to the evidence in its favor cumulated in recent 
years (see Sect. |]). The real origin of color gradients might 
actually reside in a combination of age and metallicity ef- 
fects. However, we can rather safely exclude that metal- 
licity increases in the outer parts of galaxies, because age 
variations cannot explain the color gradients when a con- 
stant metallicity is assumed (Case aA) and are just viable 
when large central ages are adopted (Case mA15). We can 
also exclude that more than half of the observed color gra- 
dients is driven by age variations, if the present day central 
ages of cluster ellipticals are < 10 Gyr. Scenarios where 
metallicity "overdrives" the color gradients by compensat- 
ing age increases in the outer regions of the galaxies are 
also possible. Approximately 0.1 dex additional decrease 
in metallicity is needed to compensate 2.5 Gyr increase in 
age every radial decade. 

In principle the use of both the H/? (more sensitive to 
age variations) and the Mg2 (more sensitive to metallicity 
variations) line indices could allow to disentangle age and 
metallicity effects (Worthey 1994). Following this idea, 
studies of the H/3 and Mg2 line index gradients in field 



(Gonzalez 1993) and in Coma cluster ellipticals (Mehlert 
et al. ^OOOa ) find on the mean that the outer regions of 
early-type galaxies are metal poorer (0.1 dex per decade in 
radius) and slightly older (0.04 dex per decade in radius) 



than the inner parts (Mehlert et al. 2000b, see also Bressan 



et al. 1996 and Kobayashi & Arimoto 1999). Our model 



eled u nder m A15. Tamura et al. (200C) and Tamura and 
Ohta ( 200C| ) reach the same conclusion from the analy- 
sis of the color gradients of ellipticals in the Hubble Deep 



ing of color gradients requires steeper metallicity gradients 
(w 0.3 dex per decade) to compensate for these kinds of 
older stellar populations in the outer regions of galaxies. 
Note, however, that the modeling of line indices is uncer- 
tain (Maraston, Greggio & Thomas 1999). Moreover, the 
possible presence of old, metal poor stellar populations 
complicates the interpretation of Balmer lines as age in- 
dicators (Maraston & Thomas 2000 ). 

It is not clear that stronger conclusions on the origin 
of color gradients of cluster ellipticals could be achieved 
looking at more distant clusters. On the one hand, even 
the resolution of HST might not be enough to allow their 
reliable determination down to the same limiting magni- 
tude at larger redshifts. On the other hand, we may likely 
fin d trac es of residual star formation (e.g., Poggianti et 
al. 199^) that make the interpretation of the evolution 



impossible, adding too many free parameters to the prob- 
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lem. Dealing with lower redshift cluster might also not be 
very informative, because the expected evolution is small. 
Improved constraints can instead be obtained by enlarging 
the sample of local and z !=s 0.4 cluster galaxies, probing 
the brighter and fainter end of the luminosity function. In 
particular, fainter ellipticals have bluer colors which might 
require too young ages and rule out the age driven origin 
of color gradients completely. 
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